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I. INTRODUCTION 



The goal of this paper is to study some analytical aspects of gauge theories in the con- 
formal window related to explicit breaking of the chiral symmetry. 1 have in mind SU{N), 
SO{N) and Sp{N) gauge theories, whose fundamental degrees of freedom are gluons and 
minimally coupled Dirac fermions (flavors) transforming under the gauge symmetry in a 
generic representation of the gauge group. 

For few flavors, chiral symmetry is expected to be spontaneously broken in such theories 
by the fermion condensate (ipip)- At large enough number of flavors, asymptotic freedom 
is destroyed. It is commonly believed that a window [conformal window) for the number 
of flavors exists in between these two regimes in which chiral symmetry is restored and 
the theory becomes asymptotically scale invariant at large distances. The existence of the 
conformal window has been conjectured by Banks and Zaks who showed that the two- 
loop running coupling constant displays an infrared (IR) fixed point if the number of flavors 
is properly chosen. The computation of Banks and Zaks can be trusted only in the large- 

Veneziano limit of an SU{N) theory. If Nf is the number of flavors of fundamental 
fermions, then the ratio Nf/N is kept fixed, and it is a continuum parameter that can be 
tuned arbitrarily close to its value at which asymptotic freedom is destroyed. The existence 
and size of the conformal window have been investigated by means of several analytical 
methods j^^S], based on educated but often uncontrolled approximations. 

Far away from the Bank-Zaks fixed point, these theories become strongly coupled and 
analytical investigation becomes inherently difficult. However scale invariance is a quite 
powerful tool. The anomalous Ward identity for the scale transformations allows to derive 
very general and interesting relationships, also when scale invariance is broken in a controlled 
fashion by dimensionful parameters in the action (like a mass term or a finite volume). Even 
though the techniques proposed in this paper are widely known and used in several areas 
of theoretical physics, their potential has not been fully exploited in the particular case of 
gauge theories in the conformal window. A recent analysis can be found for instance in {q]. 
(Broken) scale invariance can be also investigated in the Hamiltonian formalism, studying 
the action of the trace of the stress-energy tensor on physical states. This complementary 



approach has been recently exploited in 
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Before moving to the main point of this paper, it is worth reporting that many numerical 
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investigations (via lattice simulations) of the conformal window has been produced in the 
past few years. ^ The interest of the lattice community has been boosted by the proposal 
that theories inside or close to the conformal window can be suitable for building walking or 



conformal technicolor models for electroweak symmetry breaking |46|-|52||. Aside this possible 
application, the study of gauge theories in the conformal window is an interesting subject on 
its own, which will eventually lead to a better and deeper understanding of strongly-coupled 
gauge theories in general. While setting up numerical simulations and interpreting data, one 
has always to keep in mind that gauge theories in the conformal window are deeply different 
from QCD and most of the wisdom and intuition we might have from QCD does not work in 
general. Analytical works are necessary in order to develop new tools for guiding numerical 
simulations. 

For few flavors, chiral symmetry breaking is signaled by a nonvanishing value of the 
chiral condensate S = —{■ipip)/Nf. The Goldstone bosons (pions) appear as a pole in the 
longitudinal part of the current that generates the chiral symmetry (isovector axial current). 
The amplitude of the pole is related to the pion electroweak decay constant Fpg."^ If the 



quarks have a small mass m, then chiral symmetry is explicitly bro^ 



a mass Mps- The Gell-Mann-Oaks-Renner (GMOR) formula 
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ien and the pions acquire 



54j relates the quark and 



pion masses, the chiral condensate and the pion decay constant close to the chiral limit: 

mE = FlsMls . (1) 

The GMOR relation is an extremely useful tool in lattice simulations for identifying the 
chiral region, and for defining a chiral condensate which is free of additive renormalization 



with Wilson fermions 
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This work investigates how the GMOR relation generalizes in the conformal window. 



45 



^ For a recent review, see [l^]- A selection of recent results follows: (ol. [l3 
^ More precisely, if = u'y^^^d is the isovector axial current and \'n,p^) is the state of a ud pion with 
momentum and relativistic normalization 

(7r,p^|7r,g^) = 2po(27r)3^3(p - q) , 

then the pion decay constant is defined as: 

(0|A^(0)|^,p^) =Fpsp^ . 
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The starting point is the integrated Ward identity for the chiral symmetry apphed to the 
isovector pseudoscalar (PS) density P = U'j^d: 

{PHx)Pmd'x = -M . (2) 

When chiral symmetry is spontaneously broken, the pion susceptibility is divergent in the 
chiral limit. The divergence is saturated by the single-pion state propagating in the correlator 
{pole dominance) . Since the chiral condensate is finite in the chiral limit, the rhs of the Ward 
identity is also divergent in the chiral limit. The GMOR relation is derived by equating the 
divergences in both sides. 

In case of IR-conformality, I will show that the pole dominance does not hold. At the 
leading order in the chiral limit, both the pole and the continuum spectrum contribute to the 
isovector PS susceptibility. Moreover the isovector PS susceptibility is divergent in the chiral 
limit only if the anomalous dimension of the chiral condensate is in the range 1 < 7* < 2.^ 
Once properly regularized, the chiral condensate behaves in the chiral limit like: 

— {ifjilj) ^ Am 1+7* , (3) 

and the rhs of eq. (E]) is divergent. In this case a matching of the divergences of the Ward 
identity is still possible, but the GMOR relation is modified schematically as follows: 

mS = FpgMpg + continuum-spectrum contribution . (4) 

If the anomalous dimension is less than 1, the isovector PS susceptibility does not diverge 
in the chiral limit. Moreover the non-analytical term is subleading in the chiral condensate. 
For theories with a m — — m symmetry one gets: 

— (iptp) ^ acjim + Am ^+1* . (5) 

This problem can be circumvented by taking the first and second derivatives of the Ward 
identity ([2]) with respect to the mass: 



(#(2/)Pt(x)P(0))e d^xdS = 

{il)tl){z)%l)^l){y)P\x)P{Q))c d'^xd^d'^z = - 



mNf 



dm"^ 



(6) 



mNf 



(7) 
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Notice that the anomalous dimension of the chiral condensate is positive, but contrained to be less than 



2 by unitarity |56l-l58| 
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If 1/3 < 7* < 1 then the integral of the 3-point function {'ilJilj{y)P^{x)P{0))c diverges in 
the chiral hmit, therefore one can obtain a generahzed GMOR relation by matching the 
divergences of the first derivative of the Ward identity Finally if < 7^, < 1/3 also 
the 3-point function is integrable in the chiral limit, and one has to consider the second 
derivative of the Ward identity ([7]). 

It is well known that the Ward identity ([2]) contains UV-divergences. Section [TTl is de- 
voted to present the basic notations used in this paper, and to elaborate a formulation of 
the Ward identity ([2]) that is free from UV-divergences. The outcome will be the master 
equation, which will be applied to IR-conformal theories deformed with a mass term for the 
fermions. The basic tool for investigating the implications of IR-conformality is represented 
by the renormalization group (RG) equations. The scaling properties of the relevant ob- 
servables (isovector PS spectral density and Dirac eigenvalue density) are analyzed in detail 
in section lllli This is a technical but crucial step for identifying the IR-divergences in the 
Ward identity (and its derivatives). The generalized GMOR relation will be obtained for 
1 < 7* < 2 in section HVl and for < 7* < 1 in section |Vl Finally a possible application to 
lattice simulations is discussed in section IVll 

II. BASIC NOTATIONS 

A. Regulated chiral Ward identity 

As already mentioned in the introduction, the chiral Ward identity ([2]) suffers of UV 
divergences. In fact the isovector^ PS correlator Cps{x) = (P^(x)P(O)) diverges like l/|x|^ 
(up to logarithms) at short distances, which gives a quadratic divergence when integrated 
over the spacetime. In the chiral limit, it may suffer also of IR divergences. In a generic 
power-law scenario in which Cps — l/|x|° at large distances, the chiral Ward identity would 
contain an IR divergence for a < 4 (notice that in case of chiral symmetry breaking a = 2). 
The goal of this section is to derive a formulation of the chiral Ward identity in which UV 
and IR divergences are separated. This will be achieved in two steps. 

The first step consists in expanding both the PS correlator and the chiral condensate 

^ Since I am interested in this paper only in the isovector PS channel, I will omit the word isovector from 
now on. 



5 



in eigenmodes of the Dirac operator. Let Xa{x) be the eigenvector of the Dirac operator 
with eigenvalue Ua (at fixed gauge background and finite volume) and let p{u, m) be the 
eigenvalue density defined as: 

p(w, m) = lim ^y^{5{uj - uja)) . (8) 

a 

The eigenvalue density p, as all the other observables, is a function of the renormalized 
fermion mass m, the renormalized coupling and the renormalization scale /iq. However 
the dependence on g and /iq is not explicitly shown unless necessary. 

The UV divergences in the chiral Ward identity can be smoothed out by cutting the 
highest eigenvalues of the Dirac operator, by means of a regulator IZu^^iuj) = exp(— w^/w^). 
The regulated PS correlator and chiral condensate are defined respectively as: 



a,b 

Jq + 

When Uoo is sent to infinity and the UV regulator is removed, the regulated Cp:g{x,m,u 00) 
becomes the full PS correlator Cps(x,m), while the regulated S7^(m,ci;oo) diverges. The 
regulated quantities defined above have the nice property of satisfying exactly the chiral 
Ward identity: 

/ Cfsix, m, u^)d'x = ^-^^"^^^^^ . (11) 
J 

The reader should refer to appendix |X] for the detailed derivation of eq. f lTT]) . 

The second step consists in splitting and rearranging the integrals over spacetime and 
Dirac eigenvalues in the regulated chiral Ward identity: 



Hdt [ d'^c'^si^,m,u^)- r4^^LM^^= 

Jt J Jq u^ + m^ °° 

= - Tdt [ d^x C%ix, m, o;^) + r 4t^^L (^)^^ • (12) 

Jo J J Q ^ + 

This formula segregates the potential leading IR divergences in the Ihs and the potential 
UV divergences in the rhs. However since the Ihs is finite when u^o is sent to infinity, then 
the UV divergences must cancel out in the rhs as well, and the quantity 

m^{m,T,u)= lim [- [ dt[d''^CfJx,m,uj^)+[ ^fp^TZl {u)diA (13) 
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is finite. Taking the cJoo ~> oo limit of both sides of eq. (fT2|) yields: 

POO r yi 

I dt I d^x Cps{x, m)--^ = AW{m, T, to) . (14) 
Jt J 2m 

where "E^ is the chiral condensate (at nonzero mass) with a sharp UV-regulator: 



2,„/ 4^d.. (15) 



The basic feature of eq. (fT^ is the segregation of the potential leading IR divergences (in 
the m — J- limit) in the Ihs. This is easily seen by inspecting the definition of AW in 
eq. fll3p : (a) the integral of the Dirac spectral density is finite in the chiral limit since the 
region around = is not included, (b) the integral of the regularized PS correlator might 
be divergent in the chiral limit, but this possible divergence is subleading with respect to 
the one in the Ihs of eq. (|T4l) because the region around t = is not included. 

B. Spectral decomposition for the PS correlator 

At nonzero fermion mass and in absence of mechanisms that suppress particle decay 
(like finite volume or large- iV limit), any correlator is expected to have a possible pole 
contribution and a contribution from a continuous density of states with masses above some 
pair-production threshold. The pole is absent if there is no stable single-particle state in 
the considered channel. If chiral symmetry is spontaneously broken, even if a small mass 
is given to the fermions, the PS meson (pion) is the lightest particle of the spectrum, and 
hence it is stable. In the conformal window, although the PS meson is still the lightest 



isovector meson (for the Weingarten inequalities [59|]), it might not be the lightest particle 
with non-trivial isospin. In the latter case, the PS meson might be unstable. 

Having said so, the Kallen-Lehmann spectral representation for the PS correlator is 

Cpsix, m) = -^^^A(x, Mis) + / Rps{s, m) A(x, s)ds . (16) 

■^"^ JSps 

Mps is the mass of the pole, and Fps is its electroweak decay constant. The poleless case 
corresponds to a vanishing decay constant. Sps is the squared threshold energy for the 
continuum part of the spectrum. A(x, s) is the Euclidean propagator of the free scalar field. 
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Plugging the spectral decomposition into eq. ( !T^ one gets: 





Aiy(m,T,a;) . 



(17) 



I will refer to this equation as the master equation. Each term of the master equation will 
be analyzed separately in sections IIVI and |Vl I will refer to the first term in the Ihs as the 
pole-term., to the second one as the R-term, to the third one as the T,-term, and finally to 
the rhs as the reminder. 

III. RG ANALYSIS 
A. General formulae 

The renormalization procedure implies the introduction of a renormalization scale no 
and the definition of suitable renormalized parameters: the renormalized coupling constant 
g and the renormalized mass m. The choice of /io is arbitrary and has no particular physical 
meaning. In all the physical observables, a change of the scale ^ fi can be reabsorbed by 
changing the parameters of the theory g g{fi) and m — )■ rh{fi). The running parameters 
g{fi) and m(/i) are governed by the RG-equations, and can be traded for the RG-invariant 
parameters A and M. IR conformality is modeled with an IR fixed point in the coupling 
constant g{0) = The reader is referred to appendix [B] for an overview of the notations. 

Let be 0{E) an observable with the following properties: (a) it renormalizes multi- 
plicatively; (b) it is a function of a physical (i.e. RG-invariant) energy E. As shown in 
appendix [Bl one can solve the RG-equation and use the existence of an IR-fixed point to 
obtain the general form: 



where Zo(/io/A) is a renormalization factor, and 7* is the anomalous dimension of O at 
the IR-fixed point. 

The RG-analysis for the mass or the decay constant of a particle is much simpler. These 
are RG-invariant quantities, hence they do not depend on the RG-scale /io once everything 
is expressed in terms of A and M. If Mx is one of those observables, its general form is: 




(18) 




(19) 
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FIG. 1: Sketch of the running couphng for an IR-conformal UV-free theory. The transition between 
the IR-conformal region (at small fi) and the asymptotic-freedom region (at large fi) happens at 
an energy scale which is proportional to A. For instance the transition energy can defined as the 
value of at which the running coupling is half of its value at the fixed point. If A is sent to 
infinity, also p, goes to infinity, the asymptotic freedom is washed away, and the running coupling 
becomes g{fi) = g^. In this limit the theory is exactly scale invariant. 

At this point one would like to sit in the cliiral point M = and see from eq. ( fT8l) that 
the observable O exhibits the power law at small energies that is typical of IR-conformal 
theories. However this cannot be inferred unless one assumes that the limits A — oo and 
M — )■ are regular. Various combination of these limits identify different regimes of the 
theory, which I will briefly describe separately. 

• Conformal theory: M — and A — t- oo. Notice that the A — > oo is the same as 
the limit g ^ g^, at fixed renormalization scale /xq (figure [I]). Asymptotic freedom 
is washed away and the theory becomes exactly scale invariant. Assuming regularity 
of the RG- invariant functions rhx and O, one obtains that all masses vanish and the 
observable O behaves like an exact power law in the energy scale E: 



• Mass-deformed UV-conformal theory: A — !■ oo and M finite. As before, the running 
coupling is tuned at the value of the fixed point, but scale invariance is broken by the 



Mx = 0, 

(O) , (O) - 

0{E) = nl' E'^°-^* 0(0,0) . 



(20) 
(21) 
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mass: 



Mx = Mmx(O) , 



(22) 



0{E) = /.J 




(23) 



Particles stay massive and all the masses are proportional to the RG-invariant fermion 
mass M. Scale invariance is approximatively restored at small distances/high energies, 
in fact the observables O behaves like a power law asymptotically in the regime E ^ 
M. Notice that if A is set to infinity, the high energy limit is formally identical to the 
chiral limit (only the ratio E/M matters). 

• IR-conformal theory: M — and A finite: 



In this limit the theory is asymptotically free and exhibits asymptotic scale invariance 
at large distances/low energies. Masses vanish, while the observable O exhibits the 
power-law behavior described in eq. ( 1211) asymptotically at low energies -E ^ A. Notice 
that if M is set to zero, the low energy limit is formally identical to the A — )■ oo limit 
(only the ratio E/A matters). 

Summarizing, eqs. , 1^ , , , (^S) s.rid 1^ are obtained from the RG equations 
in presence of an IR fixed point, under the assumptions that the beta function and the 
anomalous dimensions are regular at the fixed point, and that the UV-conformal (A — > oo) 
and IR-conformal (M — 0) limits of the observables Mx and 0{E) are finite. I will refer to 
this set of assumptions as the hyperscaling hypothesis. 

B. Specialization to relevant observables 

I want to specialize now the general formulae obtained in the previous subsection to the 
observables we are interested in. The PS mass Mps, decay constant Fps and square root of 
the threshold energy Sps are RG-invariant quantities with dimension of a mass. Therefore 



Mx = 0, 



(24) 




(25) 



10 



the general formula (fT9|) applies: 

Mps = Mmps{M/A) , (26) 

Fps = Mfps{M/A) , (27) 

Sps = MHps{M/A) . (28) 

The relation between the fermion mass m at the scale /iq and the RG-invariant mass M 
is given by (see appendix [B|) : 

■ (-) 

From the Ward identities for the chiral symmetry, one can easily compute the anomalous 
dimensions of the PS correlator {'jcps = 27), the PS spectral density (7rps = 27), and the 
Dirac eigenvalue density (7p = 7). Specializing the general formula f lTH]l to the PS correlator 
(with E = \x\~^ and do = Q), and to the PS spectral density (with E = ^/s and do = 2) is 
straightforward: 

Cpsix) = Z2(/io/A)/iMa;|-^+2,. |^^|^|^ ^ (3q) 

Rpsis) = Z2(/io/A)/i^.^^^* Rps ^) . (31) 

Specializing the general formula ( |T8|) to the Dirac eigenvalue density is slightly more 
involved, because of the dependence on the parameter u that is not RG-invariant. Since the 
eigenvalue has the same anomalous dimension as the mass, one can use a relation similar to 
eq. (I29l) in order to define a RG-invariant eigenvalue Q: 

Finally one can take eq. ( |T8|) with E = Q and do = <i and get: 

piu) = Zipo/Ah'o^'^'' P j) ■ (33) 

IV. GMOR-LIKE RELATION FOR LARGE ANOMALOUS DIMENSIONS 

Since the PS correlator falls off like in the chiral limit, the PS susceptibility is 

IR-convergent in the same limit if < 7* < 1, while it is divergent for 1 < 7* < 2.^ In 



^ The case 7* = 1 gives rise to logarithmic divergences in the PS susceptibihty, and would require a separate 
mathematical analysis. I will just ignore this marginal case for the rest of the paper. 
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this section I will show that all the terms in the Ihs of the master equation ([T71) diverge like 
j^^2-27* ^YiQ chiral limit. Since the reminder of the master equation f|T7|) is subleading in 
the chiral limit, the divergences in the Ihs must cancel each other. This condition will lead 
to the GMOR-like relation. 



A. Chiral limit of the master equation 

Let us focus on the pole term of the master equation fll7l) . Close to the chiral limit, the 
relation f l2^ between the running mass m and the RG-invariant mass M becomes: 

m ~ Zq^i2-^'M^+"^* , (34) 

where I have introduced the shorthand notation Zq = Z{^q/ h). 

Using the general formulae fl26|) and f l27|) obtained in the previous section, the leading 
behavior of the pole-term in the chiral limit is extracted: 

^Mke---- - \zlfls^ls^V'M---^* . (35) 

Both the functions fps and rhps with no argument have to be meant evaluated in zero. 

Let us focus now on the _R-term of the master equation (fTTl) . Since Rps{s, 0) oc s^"'^* for 
small values of the squared energy s in the chiral limit, the integrand goes like In the 
range of anomalous dimensions 1 < 7* < 2 the integral is divergent. In order to single out 
the divergence, one has to plug formula ( 13T|) into the R-term of the the master equation ( |T7I) . 
and to change the integration variable accordingly to = x'^s. One gets: 

.-1/2 

/■°° //c r^ps 

/ /2p5(s,m)e-'^^— ~2ZoVo^*M2-2> / x-^+^^' Rps{x,0)dx . (36) 
Jsps ^ Jo 

Again, the function sps with no argument has to be understood as evaluated in zero. Notice 

that the integral in the rhs is finite because Rps{x, 0) is finite at x = (for hyperscaling 

hypothesis). Moreover the leading term in the mass does not depend on the UV regulator 

T^^ and on the scale A of asymptotic freedom. 



3-7* 



Consider finally the S-term of the master equation (fT7|) . Since p(a;,0) oc cji+t* for small 



1-37, 



values of the eigenvalue in the chiral limit, the integrand goes like u i+t* . In the range of 
anomalous dimensions 1 < 7=,, < 2, the exponent is smaller than —1 giving rise to a 
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divergent integral. In order to single out the divergence, one has to plug formula (13T|) and 
to change the integration variable accordingly to M = xQ, where Q is implicitly defined in 
terms of uj in eq. f l32p . After some long but straightforward manipulations, one gets: 



The integral in the rhs is convergent because p(0,0) is finite (for hyperscaling hypothesis). 
Notice that the leading term in the mass does not depend on the UV regulator u and on 
the scale A of asymptotic freedom. 

B. GMOR-like relation 

Summing up, all the terms in the Ihs of the master equation are equally IR-divergent 
(like M^"^'''*) in the chiral limit, if the anomalous dimension is in the range 1 < 7* < 2, 
while the remainder of the master equation is sub leading in the chiral limit. A GMOR-like 
relation is obtained by discarding the rhs of the master equation and dividing by the pole 
contribution: 



Eq. ( 138|) is the direct generalization of the GMOR relation to the IR-conformal theory with 
anomalous dimension in the range 1 < 7* < 2. The pole dominance is not realized in 
this case: the continuous part of the PS spectrum has the same weight as the pole in the 
regulated Ward identity f ITT]) . Notice that the continuous part of the spectrum gives always 
a positive contribution to the rhs of eq. f l38|) . 

The GMOR-like relation (!38|) does not depend on A, and therefore does not distinguish 
between the asymptotically free and UV-conformal theories. In the case of the (mass- 
deformed) UV-conformal theory, the GMOR-like relation (l38ll can be obtained much more 
easily. In fact at infinite A the PS correlator becomes UV integrable since it goes like 
for small values of x (2 < 6 — 27^, < 4). Also the chiral condensate becomes UV finite 
since the Dirac eigenvalue density goes like oji+t* for large eigenvalues (| < < 1). 
Therefore the UV cutoff can be safely removed in the Ward identity (|TT|) . and eq. f p8|) 
can be shown to be valid at every mass: 




(37) 




_-l/2 



(38) 




(39) 
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V. GMOR-LIKE RELATION FOR SMALL ANOMALOUS DIMENSIONS 



A. 1/3<7*<1 

For values of the anomalous dimension in the range < 7* < 1, the terms in the Ihs of 
the master equation f|T7j) do not contain IR-divergences. The chiral limit of the regulated 
Ward identity (fTTj) yields: 

/ C^si^, 0, u:^)d'x = r ^^K^ {u)duj . (40) 
J Jo ^ 

In fact the PS correlator decays like at large distances, which is integrable for 

1-37* 

< 7* < 1. While the integrand in the rhs vanishes like u i+t^* at small eigenvalues, which 
is also integrable for < 7^, < 1. The argument used in section IIVI for large anomalous 
dimensions does not hold anymore. However one can take derivatives of the Ward identity ([2]) 
with respect to the mass. Taking for instance the first derivative, one gets eq. (EI). I will 
show that both sides of eq. ([6]) are IR divergent if the anomalous dimension is in the range 
1/3 < 7* < 1.^ Matching the IR-divergences in the two sides of eq. (j6]) will yield the 
generalized GMOR relation in this case. Of course, in order to make everything well defined, 
one should introduce a regulated version of eq. ([6]), by taking the first derivative with respect 
the mass of the master equation in the form ( |T4l) : 

dm [2m 

In complete analogy with the analysis of the reminder of the master equation (ITTj) . if any 
IR divergence is present, the derivative dAW/dm contains only subleading divergences. 
Matching the leading divergences in the rhs of the last equation, I will obtain a generalized 
GMOR relation for 1/3 < 7* < 1. 

The first step consists in identifying the leading divergence in the mass of the two terms in 
the Ihs of eq. ( 14T|) . The RG-analysis leads for the 3-point function to the following functional 
form: 

my)PK^)Pmc,m = |,|3-..|,p°tV,|3^.. ^3 M\x\,M\yl ^) , (42) 



CO p p a 

dt / d^x / d'y (#(i/)pt(a:)P(0)),^„ - 



T 



dAW 

— (™.T,^) . (41) 



^ Also the case 7* = 1/3 gives rise to IR-divergences in both sides of eq. ([6]). However these divergences 
are logarithmic and would require a separate mathematical analysis. I will just ignore this marginal case 
in the rest of the paper. 
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where the function is finite in the chiral hmit and also in the UV-conformal (A — )■ 
oo) hmit, and it is exponentially decaying at large distances for nonvanishing mass. The 
coordinate dependence in front of the function comes from the assumption that the 
theory at A = oo is invariant under the full conformal group (and not only under scale 
transformations) . 

Integrating the 3-point function in x and y, and rescaling the integration variables with 
M, one gets the following chiral behavior: 



dt / d'^ / dS{i^^{y)P\x)Pm, 



Zl^j?^' M^-^^* / d^x / d^y 



xy 



\x 



3—7* 



\x-y 



3—7* 



(43) 



For 1/3 < 7* < 1 the integral is finite and the whole quantity diverges like M^~^'^* . Since the 
rhs of eq. (14T|) is finite in the chiral limit, also the second term of its Ihs must be divergent 
like M^~^'^* , in such a way that the two divergences cancel each other: 



/ dt I d^x.Cps(x,m) 

dm 



_d_ 

dm 



dt I d^x 



2m 



d'y {tP^{y)P^{x)P{0)), 



(44) 



At the leading order one can also replace the derivatives with the incremental ratios: 



1 

m 



dt / d^yi[Cps]sub{x,m) 



m 



2m 



0{M 



1-37. ^ 



(45) 



sub 



where I use the shorthand notation [F]sub{'n^) = F{m) — -F(O). Using the spectral decom- 
position for the PS correlator in eq. ( IT6l) . and dividing all terms in eq. ( H5l) by the pole 
contribution, one readily obtains: 



2m^ 



lim — 7, TT- — — 

m^o FpgMpg [2m 



= 1 + 

- sub fpS^PS 



— 1/2 
^PS 



X 



-3+27 



'[Rps{x,0) - RpsiO,0)]dx , (46) 



which is the generalized GMOR relation for anomalous dimensions in the range 1/3 < 7* < 1. 



B. 0<7*<l/3 



If the anomalous dimension is in the range < 7* < 1/3, then also the first derivative ( BTj) 
of the master equation does not contain IR divergences. One has to consider the second 
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derivative with respect to the mass: 



T 



m, T, u) 



2 2m dm? 

(47) 



If Ax is the minimum distance between any of the pairs of operators in the 4-point function 
in the Ihs, for large Ax the 4-point function in the chiral hmit falls off like (Ax)"^^"'"'''''*, 
which makes it non integrable in the IR for any positive anomalous dimension. Repeating 
the RG-analysis for the 4-point function one can show that the first term in the Ihs of eq. ( 147|) 
diverges like M~'^"'* in the chiral limit. Again, the rhs of eq. f H7|) is IR-sub leading and it can 
be discarded in the chiral limit. The divergence of the first term in the Ihs must be canceled 
by an equal divergence in the second term: 



dm 



dt / rf^xCp5(x,m) = / dt I d^^ I d^yd^z {ipip{y)ipip{z)P\x)P{0)), 

T J Jt 



dm^ 



2m 



0{M-^^') . (48) 



I want to replace now the derivatives with finite differences. For instance one can Taylor- 
expand the the PS susceptibility, and write the first derivative in terms of the connected 
3-point function: 

dt I d'^x.Cps{x,m) 



dm^ 



T 



At / dt 



m^ 



T 



d^^[Cps]sub{x,m)-^J^ dtJd^icJd^y{i:ij{y)P\x)P{0)),,^=o. (49) 



However the 3-point function can be constrained by the non-anomalous subgroup TiiTuNf {R 
is the representation of the fermions under gauge transformations) of the axial symmetry. 
In the conformal window the axial Z^j-^^^. is expected not to be broken along with the chiral 
symmetry. If TrNj is an integer^, then the transformation ip — i'-i^ip is a symmetry in the 
chiral limit. The 3-point function {ilJilj{y)P'^{x)P{0))c,m=o is odd under this symmetry, and 
therefore it must be zero. The previous equation reduces to: 



dm 



poo p 1 poc p 

dt j d^xCps{x,m) -^J^ dt j d^^\Cps\sub{x,m) . (50) 



^ This happens for instance for an even number of fundamental or (anti)synimetric two-index fermions, or 
for any number of flavors of adjoint (Dirac) fermions. 
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For a nonvanishing fermion mass the transformation if) — ^75^ is still a symmetry if com- 
bined with m — 7- —m. This gives for instance that the Dirac eigenvalue density is even under 
mass-sign flip, while the regulated chiral condensate is odd. This yields for the second 
derivative of S^j, / (2m) the following asymptotic equality: 



92 



dm? 



2m 



m^ 



2m 



(51) 



since the first derivative vanishes in the chiral limit. The condition ( 148|) for the divergences to 
cancel each other can be hence rewritten by replacing the derivatives with finite differences: 



— / 

m^ Jrp 



£^[Cps\sub{x,m) ~ — 



2m 



0{M- 



-47, ^ 



(52) 



sub 



which is completely equivalent to the corresponding eq. ( 145|) for anomalous dimensions in 
the range 1/3 < 7* < 1. Therefore one gets for < 7* < 1/3 the same generalized GMOR 
equation ( l46ll as for the 1/3 < 7* < 1 case. 



VI. GMOR RATIO ON THE LATTICE 



In the past few years, many lattice investigations have been devoted to the study of the 
conformal window. It is natural to ask whether the generalized GMOR relations can be 
used in lattice simulations to discriminate IR-conformality from chiral symmetry breaking. 

In case of IR-conformality, the two GMOR relations (!38|) and (146|) involve the PS spectral 
density, which is not easy to be computed on the lattice. However one can focus on the ratio 
mTja,/{FpgMpg) (that I will call GMOR ratio). I have shown in sects. HVl and IVl that the 
GMOR ratio is always greater or equal to one, and its value discriminates between the 
following three situations: 

\imm-^o ^/'".IT'* = 1 for spontaneously-broken chiral symmetry, 

\ liuim-^Q = 00 for IR-conformality with < 7* < 1, (53) 

1 < linim^o < 00 for IR-conformality with 1 < 7* < 2. 

I remind that the chiral condensate Tj^im) appearing in the GMOR ratio is regularized with 
a renormalized sharp cutoff u in the eigenvalue space, and hence it is a well-defined finite 
quantity in the continuum theory at any mass. 

Let us focus now on the lattice-discretized theory. All the dimensionless quantities defined 
in the lattice theory will be denoted with a hat. For definiteness I will consider Wilson 
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fermions.^ Since the Wilson-Dirac operator Dw is in general nondiagonalizable, it is useful to 
introduce the positive operator Hw = Dy^Dw- Let ]P(a)) be the projector on the eigenspaces 
of Hw corresponding to eigenvalues smaller than rh? + w^, being rh the quark mass defined 
via the PC AC relation. The discretized version of S^(m) is: 

S(^,m,w) = -^(tr^— , (54) 

V H\Y 

where g is the bare coupling constant, and V is the volume in lattice units. The discretized 
chiral condensate S (q, rh, u) is computable on the lattice (for an implementation of the 



projector P(iu;), see 60|). How do you take the continuum limit? One needs a dimensionful 
quantity to set the scale (for instance the Sommer radius ro)^, and the PS mass Mps to set 
the fermion mass. If fg and Mps are respectively the Sommer radius and the PS mass in 
the discretized theory, by requiring that: 

roMps{g,m) = r^Mps , (56) 

one selects a trajectory rh{g) in the parameter space. This trajectory defines how you have 
to change the bare finite mass rh while going to the continuum limit {g — t- 0) in order to 
keep the physical (renormalized) fermion mass finite. 

Before taking the continuum limit, one has to give also a similar prescription for the 
eigenvalue cutoff u. One can consider for instance the mode number per unit volume: 

C{g,m,co) = ^{tiF{co))g,rh. (57) 

The authors of [gOI showed that this is a RG-invariant quantity also with Wilson fermions. 
Once uj is renormalized, the continuum limit can be taken without extra renormalization 
factors. One can reverse the argument in order to define the renormalization of u. Choose 
a value ( for the mode number per unit volume in the continuum. By requiring that the 



The GMOR relation with Wilson fermions is already discussed in [55[ in the case of broken chiral 
symmetry. 



9 



The Sommer radius rg 61| is defined from the static potential V{r) as 



r^^(ro) = c. (55) 

If the constant c is suitably chosen, the Sommer radius is well defined and finite also for an IR-conformal 
theory. 
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dimensionless combination r^^ is fixed at any finite lattice spacing (and therefore in the 
continuum hmit): 

r',ag,m{g),u)=r^oC , (58) 

one selects a trajectory u{g). This trajectory defines how you have to change the bare 
eigenvalue cutoff u while going to the continuum limit {g — ?■ 0) in order to keep finite the 
mode number per unit volume in physical units. 

The continuum limit of the GMOR ratio is readily taken: 

j.^^ m{g)±{g,m{g),u{g)) ^ mE^cpM .^g. 

This procedure selects implicitly a value for the renormalized cutoff u that depends on the 
chosen value of C,. The relation between u and C, can be written by taking the continuum 
limit of fISS]) and writing the mode number per unit volume in terms of the eigenvalue density 
distribution: 

p{u,m)du = (. (60) 

However we are not interested in the particular value of u, since I have already shown that 
the GMOR ratio does not depend on it in the chiral limit. 

Summarizing, the GMOR ratio is a well-defined quantity and can be defined also on 
the lattice (even with a nonchiral formulation of the fermions). Its behaviour in the chiral 
limit allows to discriminate between chiral symmetry breaking, IR-conformality with small 
anomalous dimensions, and IR-conformality with large anoumalous dimensions. 

VII. CONCLUSIONS 

The well-known GMOR relation for QCD-like theories that break chiral symmetry spon- 
taneously 

hm = 1 (61) 

relates the PS mass and decay constant to the fermion mass and chiral condensate close to 
the chiral limit. It is worth noticing that all the quantities appearing in the formula above 
are well-defined: in particular the chiral condensate Sij(m) has been defined in eq. (fT5|) 
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with a renormalized sharp cutoff in the eigenvalue space. The GMOR relation comes from 
matching the IR divergences arising in the chiral limit in the two sides of the chiral Ward 
identity for the PS susceptibility I studied a possible generalization of the GMOR 
relation to the case of theories in the conformal window. Several lessons can be learned from 
the computation. 

It is possible to identify two ranges for the anomalous dimension 7^, of the chiral conden- 
sate at the IR fixed point, for which the considered gauge theories behaves quite differently: 
large anomalous dimension 1 < 7,, < 2 and small anomalous dimensions < 7* < 1. 

For a large anomalous dimension the PS susceptibility is IR-divergent in the chiral limit, 
and the (UV-regulated) chiral condensate scales like: 



One has to be very careful in taking the chiral limit of diverging integrals, because the chiral 
expansion and the integral do not commute. A clear example is eq. f p7|) : the leading term 
of the chiral condensate in the chiral limit does not depend only on the Dirac eigenvalue 
density in the chiral limit, but keeps track of the details of the mass-deformed theory. The 
correct leading behavior of the PS susceptibility and chiral condensate can be extracted only 
after studying the effects of asymptotic scale invariance on the PS correlator and the Dirac 
eigenvalue density. This is where the RG equations come into play: IR scale invariance is 
modeled as an IR fixed point for the RG flow, close to which the theory has a certain number 
of regularity properties {hypers caling hypothesis). One learns that the pole dominance in 
the PS susceptibility, which is valid for spontaneously broken chiral symmetry, does not 
hold in the conformal window. A GMOR-like relation is obtained as for chiral symmetry 
breaking by matching the IR divergences in the two sides of the chiral Ward identity, but it 
is modified by the contribution of the continuum part of the spectrum of the PS channel in 
the mass-deformed theory: 



For small anomalous dimension the PS susceptibility is finite in the chiral limit, and the 
(UV-regulated) chiral condensate scales like: 



3—7^ 

1+7* 



(62) 




--1/2 



(63) 



E, 



3 — u 



(64) 
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Considering the same ratio as in the GMOR relation, one gets: 

hm ^„ , ^„ = +00 . (65) 

However this result is completely driven by the cut-off dependent part of the chiral conden- 
sate, since the term that scales with the anomalous dimension is subleading. A GMOR-like 
relation can be obtained by explicitly subtracting the linear term from the chiral conden- 
sate. This is done by considering derivatives with respect to the mass of the chiral Ward 
identity, which involve 3- and 4-point functions. Asymptotic scale invariance is not sufficient 
for determining the IR divergences of the derivatives of the chiral Ward identity. One has 
to assume that the theory that sits on the IR-fixed point is not only scale invariant, but 
also conformal invariant (i.e. invariant also under special conformal transformations). It is 
worth to notice that the existence of an IR fixed point is related to scale invariance, but does 
not generally imply conformal invariance. However under this extra hypothesis one gets a 
GMOR-like relation for the case of small anomalous dimension: 



m 

lim 



2 r„ . . „.-l/2 



m^O FpgMpg L 



fpgrhpa Jo 



-PS „ 

^_3+27. [Rp^i^x, 0) - RpsiO, 0)]dx , (66) 



' PS"''PS 

where the sub subscript indicates that the linear term in the mass has to be subtracted by 
the chiral condensate. 

What are the possible applications of the found relations in lattice simulations? I have 
shown in sec. I VI I how the GMOR ratio ttlLq/ [Fp^M^tg) can be properly defined on the 
lattice with Wilson fermions. This quantity is interesting because it has to be 1 in the chiral 
limit if chiral symmetry is spontaneously broken, and strictly larger than 1 in the conformal 
window (in particular it has to be infinite for small anomalous dimension). Understanding 
whether the computation of the GMOR ratio on the lattice is also practically viable goes 
beyond the goals of this paper, and will be the subject of future investigations. 
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Appendix A: Regulated chiral Ward identity 



The goal of this appendix is to derive the regulated Ward identity in eq. ( ITTj) with 
the definitions in eqs. and fITU]) . The starting point is the PS correlator in Euclidean 
spacetime. Since the fermion fields enter quadratically in the action, they can be explicitly 
integrated out in the functional integral: 

Cps{x,m) = - {d-f5u{x)u-f5d{0)) = 

= (tr [iP + m)-\x, 0)75(-^ + my\0, x)-/^]) = 

= (tr [{p + my\x,0){- p + my\0,x)]) . (Al) 

In the last line, I have used the property 75 pj^ = — p. 

In a finite box with volume L'^, the massless Dirac operator at fixed gauge background has 
a discrete spectrum. In Euclidean spacetime it is antihermitean, hence it can be diagonalized 
and its eigenvalues are purely imaginary. Let Xa{x) be the generic eigenvector (spin and 
color indices are not explicitly written), corresponding to the eigenvalue iua- 

pXa = lUaXa . (A2) 

The eigenvectors can be chosen to be orthogonal and normalized: 

J Xl{x)Xh{x)d^x = 6ab . (A3) 
Both Xa and Ua depend on the gauge background. The inverse of the massive Dirac operator 



is: 



a 

Plugging this formula in the PS correlator, one gets: 

CM.,n) = lim J2^mMpM^^ . (A5) 

a,b 

The regulated PS correlator is obtained by introducing a regulator 7^^_^(a;) which cuts 
the eigenvalues above cUoo- One can choose for instance 7l^^{uj) = exp(— w^/ci;^). When the 
cutoff is removed, the full PS correlator is recovered. 



Mn^^iuJb)) . (A6) 



a,b 



lim C'^s{x,m,Uoo) = Cps{x,m) . (A7) 

OJoo— S>00 
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The regulated PS susceptibility is obtained by integrating the regulated PS correlator over 
spacetime. One can use the orthogonality of eigenvectors and the translation invariance to 
get: 

hm — dxdy > (^^ — ^ w . — — ^^a;^(wa)7^c.^(w^,)) = 



a,b 



= limi-V(^^-^^). (A8) 

a " 

The last line of the previous equation is related to the chiral condensate. In fact integrating 
out the fermion fields explicitly, using the representation ( 1A4I) and the reality of the chiral 
condensate, one can show that the chiral condensate is: 

lim — / d^x(uu{x)) = — lim — / d'^x{tT{p + m)~"^(x,x)) = 

L— >oo / L— 5>oo / 



- hm / ) . (A9) 



a " 

It is well known that this expression is UV-divergent. One can define a regulated (absolute 
value of the) chiral condensate as: 



Putting this together with the expression for the regulated PS susceptibility, one finally gets 
the regulated Ward identity: 



/ 



d'^x Cpc{x,m,Uoo) = • (All) 

m 



Appendix B: RG equations 



General formulae 



In this appendix I want to derive explicitly eq. (1181) from the RG equation: 
d d d ] 



23 



The first step for solving tlie RG equation is to introduce tlie running coupling g and 
mass m, and the renormalization factor Zq, that satisfy the following equations {RG flow): 

/^^»^(/^) = -7(^(/^))^(/^) , (B3) 

/x-^ log Zo(/i) = -lo{9{^^)) , (B4) 

gifJ'o) = 9 , rh{l^o) = m , Zoifio) = 1 . (B5) 

Using the derivative chain-rule, one can easily show that the RG equation is equivalent 
to the following one: 

^ [Zoifi)0{E, ^(/i), m(/i), ^)] = . (B6) 
This equation is readily integrated in terms of the initial conditions of the RG flow: 

Zo{fi)0{E, g{ij,), m{iJ,),n) = 0{E, g, m, /xq) • (B7) 



Before proceeding further, one needs to go back to equations of the RG flow. Eq. flB2p 
for the running coupling can be solved implicitly. One can use asymptotic freedom and the 
perturbative expansion for the /3 function: 

m = -f3og' + f3ig' + 0{g') , (B8) 

in order to write the solutions of eq. flB2p as function of a single integration constant A 
(which depends on /iq and g): 

A 1 4 f /-sM / 1 1 /?i \ 1 

-.e-^,(,)^exp|-/ (_,_,|L),,[. ,B9) 

The running coupling depends on the scale fi but also on the initial condition via the 
integration constant A. When needed, I will use the notation g{fi]A). Eq. f lB3p for the 
running mass can be integrated explicitly. However also in this case the initial condition 
can be traded for an integration constant M (RG-invariant mass) that satisfies fh{M) = M. 
The running mass can be written as: 

f /•9(m;A) ( \ ^ 
m{fi; A, M) = exp <^ - / -wAdz } M . (BIO) 
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Notice once again that the running couphng and mass do not depend exphcitly on /io once 
they are expressed in terms of the integration constants A and M . Also eq. flB4p for the 
renormahzation constant can be integrated exphcitly: 

Zo(/x; /io, A) = exp I - T^'^'^^ I . (Bll) 

[ i9(/.o;A) P\Z) J 

However a residual dependence on the renormahzation scale /io survives in this case. 

One can go back to eq. (IB7p . choose fi = E, and write explicitly the dependence on /aq, 
A, M: 

0{E, g, m, /io) = Zo{E; /io, A)0{E, -g{E; A), m{E; A, M),E). (B12) 

Being do the engineered dimension in mass of the observable O, and being Zp dimensionless, 
from dimensional analysis one gets: 

0{E, g{E; A),m{E; A, M),E) = E^oQ |^ , (B13) 
Zo{E■^^,.^) = Zo(^,^^ . (B14) 



^ A ' A^ 

Plugging these last two equations into eq. (]B12p . one finahy gets: 

0(E,,,m,/io) = Zo(^,|)i^^°o(f,|) . (B15) 

2. Mass-deformed IR-conformal theory 

IR scale invariance is recovered in the chiral limit if the beta function has a zero in g = g^, 
and the coupling constant runs into it in the IR limit g{fi ^ 0) = g^. I will assume regularity 
of the beta function and the anomalous dimensions around g = g* (hyperscaling hypothesis): 

m = ^{9-9*) + 0{{g-g.f) , (B16) 
l{9) = l* + 0{g-g,) , (B17) 
lo{9)=lf^ + 0{g-g,) . (B18) 

The goal of this subsection is to specialize the formulae of the previous one to this particular 
case. 

The running coupling can be expanded at small renormahzation scale /t <^ A: 

^(/i)^^/.-A(£)''* . (B19) 
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Notice that /3* > since (7* is an IR fixed point, and A > since we are considering 
asymptotically-free solutions with g{oo) = 0. 
It is useful to introduce the following function: 



Zoifi/A) 



JO) 



JO) 



9* - gjf^/A) 

A 



exp 



io{z) ^ r* 



dz 



(B20) 



which is normalized in such a way that ^o(O) = 1, and in terms of which the renormalization 
factor is written as: 



(O) 



One can plug the Zq function into eq. flB15l) . and define O = 0/Zo{E/A) in order to 



get: 



0{E,g,m,^io) = Zo (^) f^rE'^o-.^'o 



(o) . /M E 



(B22) 



Replacing 70 with 7 in eqs. (IBllI) . (IB 141) . ( 1B20I) and ( ]B21I) . one can define the quantities 
Z, Z, Z in analogy with Zq, Zq, Zq- The running mass defined in eq. ( 1B10I) can be written 
as: 



(B23) 



which displays the typical power law of IR-conformal theories. Notice that Z is the renor- 
malization factor of the chiral condensate (therefore the mass renormalizes with Z~^). 
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